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The response of an extended periodic system to a homogeneous field (of wave- vector q = 0) cannot 
be obtained from a q = time-dependent density functional theory (TDDFT) calculation, because 
the Runge-Gross theorem does not apply. Time-dependent citrreni-density functional theory is 
needed and demonstrates that one key ingredient missing from TDDFT is the macroscopic current. 
In the low- frequency limit, in certain cases, density polarization functional theory is recovered and 
a formally exact expression for the polarization functional is given. 
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Density functional theory0,0 is a standard approach 
for calculating ground-state properties of solids Q and 
molecules Time-dependent density functional theory 
(TDDFT) is an extension of the ground-state formalism 
based on the Runge-Gross theorem this establishes a 
one-to-one correspondence between time-dependent den- 
sities and time-dependent one-body potentials. When 
a time-dependent electric field is applied to a system, 
this formalism provides a route to its optical response @. 
The response equations of TDDFT have been encoded 
in standard quantum chemical packages 0, and results 
for molecules are appearing (see Ref. [8j for many exam- 
ples). As in the ground-state case, the accuracy depends 
on the quality of the approximate functional used. 

There is great interest in applying the same tech- 
nique to extended systems. While these can be treated 
well within existing wavefunction technology, using, e.g., 
the GW approximation and then solving the Bethe- 
Salpeter equation for the optical response^, the al- 
lure of a TDDFT approach is its far lower computa- 
tional cost. Calculations already show that excitonic ef- 
fects appear to be treatable by going beyond the usual 
local and semi-local approximations of standard DFT 
calculations ^3 El ■ 

There is also a version of the time-dependent theory, 
called time-dependent current-density functional theory 
(TDCDFT), that uses the current-density as the ba- 
sic variable: As the choice of variable (charge den- 
sity versus current density) appears a matter of con- 
venience, TDCDFT and TDDFT appear to be equiv- 
alent (for non-magnetic systems). The time-dependent 
exchange-correlation potential has been argued to be 
more amenable to local and semilocal approximation in 
terms of the current-density 0] an d this framework has 
been used in recent response calculations of solids 
and conjugated polymers [lfij . Initial work towards a 
matrix formulation of the current-density response equa- 
tions has been presented in Ref. ^(| . 

In this paper, we demonstrate a difference in princi- 
ple between the two approaches when applied to bulk 
solids. The basic theorems of DFT, ground-state or time- 
dependent, are proven for finite electronic systems (i.e. 
systems with a boundary). We consider the response 
of periodic systems (such as the bulk of an insulator or 



metal) to time-varying electric fields which have a spa- 
tially uniform component. We show that TDDFT fails 
in this case: there is no one-to-one correspondence be- 
tween the time-varying periodic density and the applied 
potential. Instead, TDCDFT is needed for a complete de- 
scription. The time-dependent potential and all response 
properties are not functionals of the time-dependent pe- 
riodic density, but rather are functionals of the time- 
dependent periodic current density. In the limit of low 
frequencies our analysis recovers the well-known GGG 
theorem |l7j . if microscopic transverse currents can be 
neglected; TDCDFT then recovers (static) density po- 
larization functional theory [ItI Hsj . 

Modern solid-state calculations model extended peri- 
odic systems and extract bulk properties by using pe- 
riodic boundary conditions. Our first point is that the 
Runge-Gross (RG) theorem does not apply when a ho- 
mogeneous electric field is applied to a periodic system. 
RG states that, given an initial state, there is a one- 
to-one correspondence between time-dependent densities 
and time-dependent scalar potentials for a given inter- 
action and statistics. The first step of the RG proof 
establishes a one-to-one correspondence between poten- 
tials and currents ,20). In the second step, continuity 
is then used to relate currents to densities, and a one- 
to-one mapping between densities and potentials results 
provided a certain surface integral involving the density 
and the potential vanishes. For finite physical systems 
the condition for requiring this surface term to be zero 
can be given rigorously for systems in which the density 
vanishes at the surface j^. For a periodic system, one 
might try to choose a surface around which the density 
and potential are periodic but for a uniform field the lin- 
earity of the potential prevents this, and TDDFT does 
not apply. 

Another way to see this is in modelling the periodic 
system in an electric field by a large ring of length L of 
the material, through which is threaded a time-varying 
solenoidal magnetic field [Tj| . This field produces a uni- 
form vector potential on the system, A(t), that corre- 
sponds to a macroscopic electric field along the ring. The 
beauty of this appoach is that the Hamiltonian remains 
spatially periodic, although it becomes time-dependent. 
TDDFT was however derived with only scalar potentials 
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in mind, and does not consider such transverse vector po- 
tentials; such uniform electric fields cannot be generated 
by a charge distribution. The first part of RG still holds 
in this case however, showing the potential is a unique 
functional of the current-density. 

Thus the density in the interior of any system is insuffi- 
cient information to deduce the external electric field, but 
the current-density is sufficient. This is a time-dependent 
generalization of the orig inalGGGtheorem[l7|. We shall 
come back to the static case shortly. 

A simple example demonstrating the non-uniqueness 
of the potential-density mapping is a noninteracting free 
electron gas on a ring, subjected to a constant, uniform 
electric field, £ , turned on at time t = 0. Representing 
the field by a vector potential A = —c£t, each orbital in 
the Slater determinant state satisfies the time-dependent 
Schrodinger equation 

(p-£t) 2 (t> m /2 = i<pm (1) 

(We use atomic units throughout this paper). If the 
electrons are initially in an eigenstate, <p m (x,0) = 
e l2 ^ mx / L /t/L, where x is the coordinate around the ring, 
and the conjugate momentum is k rn = 2-7rm/L, with m 
an integer between and L, different for each orbital, 
then the solution at time t is readily found to be 

<t> m {x>t) = e -i(C*/2-fc m ^ 2 /2+S 2 t 3 /6) e ^mx/L / ^/ I> (2) 

Since the electric field only affects the phase of this or- 
bital, its density, and that of the noninteracting gas, 
n(x,t) = J2m=i \4>m{x, t)\ 2 , remains spatially uniform 
forever. In particular, two different electric fields give rise 
to exactly the same time- (in) dependent density. Thus 
the external potential is not uniquely determined by the 
density here. (This argument holds for any number N of 
electrons). 

The first part of RG remains valid 2(| and applies 
to arbitrary vector potentials, not just those describing 
an electric field. Choosing a gauge in which the scalar 
potential vanishes, one can show there is a one-to-one 
correspondence between A(rt) and j(ri) for a given ini- 
tial state [52, 13 an d this provides the formal basis 
for TDCDFT. In this gauge all electric fields are rep- 
resented by vector potentials but the one-to-one corre- 
spondence is of course gauge-independent. In our sim- 
ple example, the physical current density is given by 
= El km + N£t/L . In two different electric fields, 
two different currents arise. More generally, for a periodic 
potential on a (1-dimensional) ring, when a uniform elec- 
tric field is turned on, the density, current, etc., remain 
periodic, and each can be written as ^2 ja exp(iGx), 
where G — 2nir/a, and n is an integer. All components 
of the current-density at G ^ are determined by the 
periodic density by the continuity equation: 

j G (w) = wncM/G, (G^O). (3) 

When a uniform electric field is present, the G = com- 
ponent (the macroscopic current) is undetermined by the 
time-dependent periodic density. 



When the wavelength of the external field is finite, RG 
does apply to the periodic system; the macroscopic cur- 
rent is a functional of the periodic density (although can- 
not be determined by Eq. 10)- Examining the second 
step in the RG proof in one-dimension, one may choose 
the surface to be two points separated by an integral mul- 
tiple of the lattice constant and the wavelength of the ex- 
ternal field: the Hamiltonian is periodic around this sur- 
face, and the integrand of the surface term is the same at 
the two points. Thus the surface term vanishes and the 
RG theorem holds. This is consistent with the ground- 
state case 01 where the macroscopic polarization is only 
an independent variable when the external field has a ho- 
mogeneous component. In our simple example, consider 
applying a perturbative electric field £ cos(qx) around 
the ring, where a finite number of wavelengths fits into 
the ring so q is an integer multiple of 2n/L. The density 
picks up a spatial modulation, which, in the limit that 
q — > 0, becomes, to first order in £ , 

n(xt) -> N (l + £q ^ (x + k m t) 3 /3 ) /L. (4) 

V m=l / 

Different fields yield different densities except when q = 
0, consistent with the 1-1 mapping between densities and 
potentials at finite wavelengths. One can then imagine 
attempting to find the uniform-field value of certain re- 
sponse functions |25| by taking the q — > limit of a series 
of finite-g TDDFT calculations, although the larger su- 
percell required might render this procedure impractical 
(although, see discussion of Ref. [l(| below). The same 
issue arises in the ground-state case, where it was shown 
in Ref. [2(| how careful use of a sawtooth potential on 
a supercell can resolve this problem. The approxima- 
tion for the density-functional must be ultra-nonlocal in 
space in order to capture the exchange-correlation fields 
generated from the charge distribution modulated by the 
long wavelength. The first excitonic peak in the dielec- 
tric function of silicon (missing when using any local 
or semi-local approximation in TDDFT) was captured 
in Ref. ^lj in a TDDFT calculation at finite wavevec- 
tor, using an ultranonlocal kernel / xc ~ o/q 2 however 
with a being empirically determined, and with apply- 
ing a GW shift to the spectrum. In Ref. [ljj an exact 
exchange calculation was performed, reproducing the ex- 
citonic peak of silicon, without having to battle through 
a monstrously large numerical computation. The authors 
were able to analytically extract the diverging ultranon- 
local factor 1/q 2 from the exact-exchange kernel and so 
needed to numerically evaluate only the remaining ma- 
trix elements, well-behaved as a function of q: because 
the long-range behaviour had been already factored out, 
the numerical computation converged rapidly with the 
number of Appoints. 

However, like density-polarization functional theory in 
the ground-state case, TDCDFT allows one to obtain 
the homogeneous field result directly from a single q = 
calculation; this cannot be done in TDDFT. 
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Having established that the current is needed when 
calculating bulk response properties in solids in uniform 
fields, we next review how such a calculation is performed 
for any electronic system. The one-to-one correspon- 
dence between A(ri) and j(ri) can be used to establish 
a set of KS equations in which noninteracting electrons 
move in a KS vector potential A s (r£) and reproduce 
the exact current-density j(rt). Assume the exact KS 
ground-state has been found and is nondegenerate. The 
TDCDFT KS equations are 



p + -A s [j](rt) 



+ u s (r) } (f>i(rt) = i 



at 



(5) 



where v s (r) is the periodic ground-state KS potential 
and all time-dependence is in the KS vector potential 
A s [j](rt). The orbitals begin as the occupied ground- 
state KS orbitals. The KS vector potential is defined to 
produce the exact physical current density 

j(rt) = 3Y>*(rt)V&(rt) - -n(vt)A s (vt) . (6) 

IOCC 

Departing from convention, we write 

As(ri) = A ext (rt) + A EM (rt) + Axo(rt). (7) 

Here Aem (rt) = A H (rt) + At (rt) is the full electromag- 
netic potential, satisfying Maxwell's equation: 



1 d 2 



An . 



V r A em- V- (V- A EM ) = — r j. (8) 



The longitudinal component, A H , is the vector equivalent 
of the Hartree potential, while the transverse component, 
At, arises from the transverse component of the current: 



A T ( rt ) = l/dy ixM-lr-r'l/c) 



(9) 



In previous formulations of TDCDFT (e.g., Ref. p^. 
this term does not appear explicitly. Since At is a nonlo- 
cal classical electromagnetic contribution of any moving 
charge density, it is not an exchange-correlation effect, 
and should be included exactly in any time-dependent 
calculation. Unlike the Hartree contribution, it is not 
adiabatic, i.e., it depends on the retarded current, not the 
instantaneous density. Consider the one-dimensional ex- 
ample above cast in three dimensions: we fatten the ring 
in the radial and cylindrical directions, imposing either 
hard wall or periodic boundary conditions in these two 
directions. The current flows only azimuthally and uni- 
formly on the ring and so is purely transverse in the sense 
that it has nonzero curl and zero divergence. At the same 
time, we call it macroscopic, since it is uniform along 
the ring. The classical response appears purely from At 
rather than from Hartree. We have checked that the VK 
approximation for the exchange-correlation vector poten- 
tial |l2( is unaffected by the addition of At- (We note 



that the addition of At to static current-density func- 
tional theory has been discussed for example in |27j.1 

We pause to make a connection with TDDFT. The 
TD density is given exactly in TDCDFT, via continu- 
ity n(rt) = -V • j(ri) but a TDDFT calculation is only 
guaranteed to reproduce the longitudinal component of 
the current. If there does exist a TDDFT KS poten- 
tial v s (rt) that reproduces the interacting current as 
well as the density 0, then the TDCDFT potential is 
A s (ri) = c J* dt'S7v 3 (rt'). But there exist special geome- 
tries, like the ring, to which TDDFT cannot be applied 
at all when macroscopic fields are present. 

Next we consider the special case of linear response 
to a uniform electric field, A ext (rt) = c£ (cu) exp(— iujt + 
T)t)/iu), where t] = 0+ turns the field on adiabatically 
from t = — oo to 0. If the current does not grow indefi- 
nitely (so that Fourier transforms exist), all additions to 
the external potential in Eq. (J7J can be written to first 
order in the perturbation: 



A T (rw) = - / d d r' 



(iuj) 

1 f ^„/^(Ml r -^l/ c ))jT(r'w) 



A xc (rw) = / dV / xc [n ](rr'w) ■ j(r'u) 



(10) 



where / xc is a nonlocal exchange-correlation tensor func- 
tional of the ground-state density, analogous to the scalar 
exchange-correlation kernel of TDDFT. 

We return now to solids in time- varying electric fields, 
and consider the macroscopic response, using the ring 
geometry as described earlier. Let G be the reciprocal 
lattice vector. The density change when the external field 
is turned on 

5n(rt) = <5?ig(w) exp(iG • r — iut) (11) 

G#0 

has no macroscopic (i.e., G = 0) component due to 
charge conservation. This implies that the Hartree re- 
sponse remains always periodic and has no macroscopic 
component. Performing a spatial Fourier transform on 
the transverse potential yields 



k-TG 



Aire 



c 2 G- 



Jtg(w)- 



(12) 



We distinguish the microscopic (G ^ 0) transverse cur- 
rent from the macroscopic (G = 0) current, which is a 
spatially uniform current travelling around the ring. As 
u> — > 0, the microscopic contribution vanishes relative to 
Hartree, because of the factor of uj 2 in Eq. (|10fl . The 
macroscopic component, however, does not, yielding 



4nc 



JG=oM 



= E7*cogMM- jgH- (13) 



G 
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Thus, from the ring perspective, the origin of the elec- 
tromagnetic response (traditionally considered a Hartree 
effect) is the transverse potential generated by the ring 
current. In the limit when L — > oo, the distinction be- 
tween transverse and longitudinal breaks down, but for 
any real calculation with discrete fc-points in the Brillouin 
zone, this distinction is important. For any L, A™ ac (w), 
being the classical macroscopic response, cannot be ne- 
glected in considering solids in electric fields. 

For insulators and metals at finite frequencies, after 
the adiabatic turn-on of the field, the system settles into 
a steady state, thanks to the damping factor -q. Because 
we are using TDCDFT, the macroscopic current in this 
KS calculation is the physical one at all times. Thus the 
macroscopic polarization 

Pmac/i\ f ill • T»mac/ \ JG— o(^) 
(t) = - / at Jg=o(£ ), P M = : 

(14) 

equals the physical one. Multiplying the vector potential 
by ioj/c yields the KS electric field: 



E xc°M - — /xcooHP mac M + - J2 

c c — ' 

G#0 



7xcOGM-{^M+j T GM}(15) 

where we have used continuity. If the microscopic trans- 
verse currents vanish or can be neglected, then jc can 
be written solely in terms of no, meaning that the entire 
corrections to the external field can be given in terms 
of the macroscopic polarization and the periodic density 
change. Now consider the final state of the system af- 
ter the field has been adiabatically switched on. Then in 
the limit that uj — > 0, this second point is precisely the 
GGG assertion made in density polarization functional 
theory: that the exchange-correlation field depends on 
both the periodic density and the macroscopic polariza- 
tion [UBS. That result arises here ou t of TDCDFT, 
even at finite frequencies, and Eq. (|15fl then provides an 
exact expression for the polarization functional. However 
if microscopic transverse currents cannot be neglected, 
then for the GGG result to hold, these currents them- 
selves would need to be functional of the periodic density 
and macroscopic polarization; this remains to be investi- 
gated. 

During the turning-on of the electric field, in the oj — > 
limit the system remains always in its adiabatic ground- 
state. Thus minimization of the energy (including the 
exchange-correlation electric field terms) over periodic 
functions [l8L 12^. l30| , within the Berry-phase formalism 
[HSU, wil1 y ield the same result as the TDCDFT cal- 
culation in the low- frequency limit. 

This work also shows that polarization is an infinite 
memory effect within TDCDFT jsj, i.e., after the elec- 
tric field has reached a finite value, the system 'remem- 
bers' forever the current that flowed in turning it on. (It 



does not depend on the procedure in reaching the steady- 
state, rather it depends just on the time-integrated cur- 
rent). Similarly, in pumping a finite system from the 
ground-state to a given excited state, the current pro- 
vides a natural way for the system to remember indefi- 
nitely which state it is in. Thus TDCDFT may provide 
a natural solution to some of the paradoxes generated by 
a pure time-dependent DFT|3l|. 

Lastly, we discuss existing approximations and calcula- 
tions. Just like the Hartree term, any density functional 
approximation, e.g. ALDA or AGGA, misses entirely 
the macroscopic contributions discussed here[33. Even 
if we regard the bulk insulator as a large but finite slab, 
so that TDDFT does in principle apply, we would need 
an ultra-nonlocal functional of the density in order to 
capture the contribution to the XC field from the po- 
larization charge density on the surfaces. In the static 
limit, this is (often) essentially the GGG effect discussed 
above, where in addition to the bulk periodic density, the 
functional depends on the bulk macroscopic polarization. 
For example, for a large rectangular slab with translation 
symmetry, the exchange-correlation electric field is 



E xc (q,w) = -iq/ xc (q,w)n (1) (q,w) 



(16) 



where the TDDFT kernel in the long wavelength limit 
may be expressed in terms of the component of the 
TDCDFT kernel parallel to the field: lim q ^ / xc (w) = 

— w 2 / xc00 || (uj)/(cq 2 ). The ultranonlocal 1/q 2 dependence 

multiplies the response density, n^: this has a lattice- 
periodic part in the bulk as well as a part proportional to 
the polarization charge density on the surface. The fact 
that one needs to look only at a unit cell in the bulk and 
integrate up the current that has flowed through it to ob- 
tain P mac (i), shows how local approximations in terms 
of the current capture the essential physics, while local 
density approximations cannot. Since orbitals depend 
ultranonlocally on space, orbital functionals in TDDFT 
32] may capture polarization effects. However, since po- 
larization can be obtained from the current-density, these 
more complex functionals are not necessary to describe 
these effects if the current-density is used as the basic 
variable. 

TDCDFT calculations for the optical response of solids 
have already been reported [13|>[l4[ that explicitly include 
the macroscopic transverse component of Aem on top 
of the VK approximation. In adding the VK correction 
to their calculation^^ > the microscopic contributions in 
Eq. (Unj are not included (Eq. (16) of Ref. 14]). This 
may be related to the need to divide the beyond- ALDA 
contribution by a factor of 2.5. Similarly, the simple ap- 
proximation in density polarization functional theory of 
Ref. [35|], in which E xc = <yP mac ; where 7 is a constant, 
does not include such contributions. 

We conclude by noting that while the present work has 
focussed on the special case of the interior of a bulk in- 
sulator or metal, the underlying logic was motivated by 
the need to construct an approximately local theory of 
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exchange-correlation for any system in an electric field. 
Two examples make this clear. In tunnelling through 
a quantum wire, present calculations use ground-state 
DFT KS orbitals as their starting point. This leads to 
resonances at the positions of bare KS orbitals. For fi- 
nite systems without currents, regular DFT tells us there 
are significant exchange-correlation corrections. The only 
way to calculate such corrections for a quantum wire is 
using TDCDFT, in order to handle currents. In another 
area, electron and energy transfer in biological molecules, 



attempts are being made to estimate matrix elements 
in a TDDFT calculation|36|. Such calculations, using 
adiabatic local and gradient-corrected approximations, 
clearly miss any contributions from macroscopic currents. 
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